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ABSTPuACT 



The designer of stiffened cylindrical pressure vessels 
is often faced v;ith the problem of comparing structures 
of wiedly varying geometries. To do this, he makes use 
of a number of dimensionless parameters v/hose function is 
to reduce the number of variables in the design problem. 

The object of this thesis was to investigate the 
effect of the dimensionless parameter, X, on stiffened 
cylindrical pressure vessel design. In particular, the X 
region where collapse of vessels due to yield-initiated 
buckling between the frames v/as identified and delimited. 
This region was used as the first iteration in a design 
procedure aimed at producing a vessel of maximum structural 
efficiency. Design curves which may be used to develop 
such a vessel are provided for shell thicknesses up to 
h/Dj^ = .01. These curves are based on both maximum 
principle stress and Hencky-von Mises stress criteria and 
provision is made in a computer program for use of other 
methods. A modified Foppl Formula is used for frame 
evaluations . 

The formulas used range from conservative to liberal 
in their approximations of collapse pressures. However, 
use of even the most conservative formulas results in 
structural efficiencies which compare favorably with 
present-day practices. 

Thesis Supervisor: John Harvey Evans 

Title: Professor of Naval Architecture 
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NOMENCLATURE 



A Cross sectional area of frame, in^ 

Ab Cross sectional area of frame plus plating of 

width b, in^ 

A]^,A 2 ,A 2 Plasticity coefficients, see Reference (10) 
b Width of frame at shell, in 

c Distance from neutral axis of bending to extreme 
outside fiber, in _ 

D Outside diameter of shell, in 

d Depth of rectangular frame, in 

Djj Diameter to neutral axis of frame plus plating of 
width b, in 

Dm Diameter to neutral axis of shell alone, in 
E Modulus of elasticity, Ibs/in^ 

e Eccentricity, in 

Eg Secant modulus of elasticity, Ibs/in^ 
f Width of flange of T frame, in 
h Shell thickness, in 

Moment of inertia, frame plus plating of width b, 
in'* 

Ig Moment of inertia, frame plus plating of effective 
width, see Reference (10), in'* 

If Moment of inertia, frame alone, in'* 

Moment of inertia, frame plus plating of width L, in'* 

L Unsupported span of plating, in 

Lj^ Bulkhead spacing, in 

Lf Frame spacing, in 

Unsupported span of plating taken as frame sapcing, in 



n Nuraber of complete waves in buckled configuration - 
circumferentially 

Pc Collapse pressure, Ibs/in^ 

P]^ Actual collapse pressure desired, Ibs/in^ 

Pcl/Pc2/Pc3 Temporary values of p^ for iteration 

collapse load per unit axial length of a ring frame, 
lbs/in=^ 

R Outside shell radius, in 
Rj Radius to neutral axis of "frame, in 

Rjj, Radius to neutral axis of shell, in 

R]^ Radius to neutral axis of frame plus shell of width 

L, in 

V Volume displaced by cylinder per unit length, in^/ft 
Wt Weight per unit length, Ibs/ft 

A Weight of s^lt water displaced by cylinder per unit 
length, Ibs/ft 

ri Efficiency - defined in Section III 
X Slenderness Ratio - defined in Section I 

y Poissons ratio 

0p Axial stress at mid-thickness, Ibs/in^ 

0 IJI Tangential (circumferential) stress at mid-thickness, 

Ibs/in^ 

0y Material yield stress, Ibs/in^ 

\p Pressure factor - defined in Section I 
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I . INTRODUCTION 

Stiffened cylindrical pressure vessels under external 
hydrostatic pressure fail by one or any combination of 
three collapse modes. While both theoretical and 
empirical expressions have been developed to determine 
collapse pressures for all modes, it is up to the designer 
to determine what structural propottions will result in 
the most efficient use of material strength and stiffness 
properties for a required collapse pressure. Two dimen- 
sionless parameters which have been widely used to define 
different geometries and their capabilities are and 
X. li) is called pressure factor and it is a ratio of a 
structure's actual collapse pressure to the collapse 
pressure of an unstiffened tube with the same h/Dj^ ratio. 



X is called slenderness ratio. It is analogous to the 
slenderness ratio, 1/r, of column theory, however, it 
includes terms which allow for different material properties. 



X 




Despite t.ieir fairly wide use by designers in 
comparisons of different structural geometries, this 
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researcher is unaware that any attempt has been made 
to develop a relationship between these parameters which 
would indicate their effect on structural efficiency. 

Also, it has been stated that these parameters are 
sometimes misleading (7) , however, no definite limits 
on their usefulhess have been established in the litera- 
ture, and only a few other dimensionless parameters have 
been developed to be used instead. A need exists therefore, 
to define the limits of usefulness of both \p and X and 
to develop within those limits, relationships between 
them and structural efficiency which may be readily 
applied in the early design stages of a stiffened cylindrical 
shell . 



In this report, two failure criteria, maximum circum- 
ferential stress and Hencky-von Mises, have been used 
to develop relationships between X, \pr and efficiency. 

In the course of the development, certain trends of the 
parameters have been noted and in the concluding remarks, 
useful limits of X and \p are indicated. In particular, 
the "low lambda" region has been defined and an intro- 
duction made to the concept of maximizing structural 
efficiency within that region. Design curves are presented 
to enable optimization within the low lambda region of 
structures v:i^h collapse pressures up to 2500 psi and examples 
of the process are given. Special note should be taken 
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II. DEVELOPMENT OF DESIGN CRITERIA 

Historically, the largest and most critically 
designed stiffened cylindrical pressure vessels have 
been intended for use as submersibles . Such vessels, 
under external hydrostatic pressure, have been known to 
collapse in one, or a combination of three modes. Since 
the mechanism which precipitates failure in each mode 
is different, the designer must take all three modes 
into consideration when designing for a given collapse 
pressure. Hov;ever, he has a choice in the design method 
he may use when developing a highly efficient structure. 

He may design toward one mode of collapse, seeking to eli- 
minate the possibility of collapse in others v;hile changing 
structural proportions to obtain a high efficiency. His 
alternative is to simply design tov;ard the most efficient 
structure, not worrying about the mode of collapse as long 
as the required collapse pressure is obtained. 

An example of the first method may be found in the 
design of the relatively shallow-diving submarines of 
World War II. These vessels were designed to collapse by 
asymmetric buckling between frames. It was felt that 
imminent collapse v;ould be signalled by the appearance 
of lobes girdling the hull and the crew would be given a 
chance to surface before total collapse occurred (22) . 
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This explanation indicates that for the volume-limited 
submarines of that era, efficiency criteria could be 
established to a large degree by safety considerations. 

For deeper diving, weight-limited submarines, 
the first method proved inadequate and the second method 
v;as refined in a number of ways to allow development of 
highly efficient vessels. A refinement which will now 
be discussed is one which assumed that ip is limited to 
1.0. This assumption used the maximum circumferential 
stress criterion and minimized weight by using the 
smallest possible number of frames. It resulted in a 
structure proportioned such that X = .8. In order to 
visualize this approach, figure 1 has been developed. 

Figure 1 is a plot of ip versus d/h (frame depth/ 
shell thickness) for four X values. As implied by 
the abscissa scale which indicates relative frame strength, 
the frames are rectangular and only their depth is varied 
to vary strength. The basic geometries used in this fig- 
ure are taken from a report by Trilling (19) and one of 
his curves (curve 2), as v;ell as his experimental data 
(crosses) are plotted. The other three curves plotted 
on figure 1 represent the effect of changing X. In all 
three cases,, ■‘he distance between frames was the only 
parameter varied. Values of >p were obtained by using 
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von Sanden and Gunther's formula #92A to predict collapse 
pressure . 

It will be noted that curve 3 of this figure lies 
almost entirely along the line, ip = 1.0. The X value 
of this curve is .814, and it is this relationship which 
forms the basis of the design method cited above. Just 
as the earlier method of designing for a particular 
collapse mode proved inadequate, this method has also 
been unable to develop the efficient structures necessary 
to resist high pressures. The remainder of this section 
will be devoted to developing a design method which will 
result in highly efficient structures. To begin this 
development, a further look at figure 1 is necessary. 

On the basis of curve 2 and the experimental points, 
one of Trilling's conclusions was that "except for very 
short thick tubes, . . . according to formula #92A the 

stronger the frame, the weaker the vessel." This was 
explained by noting that very strong frames form "hard 
spots" v;hen, under external pressure, they contract less 
than the shell. The hard spots, in turn, cause stress 
concentrations in the otherwise uniformly contracting 
shell and lower collapse pressures. 



Curve 4 indicates, however, that this conclusion is 
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not necessarily true, even for thin shells and large frames. 
The only difference between the model of curve 4 and 
Trilling's models is the decreased frame spacing which 
results in different A values; .898 for Tilling 's models 
(curve 2) and .683 for curve 4. According to curve 4, 
shell strength actually increases with frame size when the 
frames are relatively close together. The seeming contra- 
diction of Trilling may be explained by an extension of 
his logic. If the frames causing stress concentrations 
are closer together (low lambda) the overall stress 
pattern will be uniform enough that premature yielding 
and buckling will not take place. 

The obvious lower limit of A is obtained when the 
frames are in contact (A-^0) (24) thereby creating a thicker 

shell (rectangular frames) or a webbed sandwich-type 
construction (T-frames) . Both of these geometries result 
in putting the frame-caused stress concentrations in such 
close proximity that the overall stress pattern approaches 
uniformity. This uniformity of stress pattern implies that 
the yield strength of most of the material in the structure 
will be utilized before collapse occurs, which, in turn, 
implies that a structure of very high efficiency is possible 
in this A area. 



With a lov;er limit of A established at 0, an upper 
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limit is necessary to define the range of geometries 
where the effect of increasing frame strength is not 
detrimental to collapse pressure. This limit will be 
referred to as For structures with X values above 

Xjfj, increasing the frame strength v;eakens the vessel. 

For X values below X^^, increasing the frame strength 
strenthens the vessel. An initial value of this limit 
can be found by direct use of VJindenburg ' s approximation 
of von Mises' formula for collapse between frames due 
to elastic shell instability. The development of this 
formula is outlined in Appendix A. In its simplest 
form, this formula is: 



'P = 1-3 

X2 



Recalling that xlf = 



^~D^y 



, Windenburg ' s formula 



determines instability collapse in terms of the hoop stress 
relation, i.e. , 

« h o,, 1.3 

Pc = -57- 

therfore, v;hen 1.3 = 1.0, failure is due primarily 

X^ 

to yielding of the material rather than plastic instability. 
At this point, X = 1.14. If X 1.14, then instability 
collapse is indicated since 1.0. In this situation. 
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strengthening the frames cannot delay collapse. It may 
be tentatively assumed that when X 1.14, strengthening 
the frames increases overall strength. If this is so, 
then an initial approximation of X^ is 1.14. It is important 
to note that this is a very approximate method for 
determining X^j^. A glance of Trilling's results on figure 1 
will indicate that collapse at \p values below 1.0 has occurred 
in structures with X values as "low as .89. Also, with 
this method, Xj^ is independent of any individual structure 
parameters. Since this is questionable, another means 
of determining Xjj^ is necessary. 

Returning to curve 3 on figure 1, it will be recalled 
that for the X value used in that curve, \p = 1.0 for all 
frame sizes. As long as the frames are of adequate 
strength or stiffness to carry a given percentage of the 
shell load, their size has no bearing on the strength of 
the cylinder. Figure 2 elaborates on this point. The 
curves of figure 2 are collapse pressures, determined 
by w92A, plotted as against values of X. Each curve 
represents one size of frame and all structures have the 
same h/D' ratio, so that only the distance between frames 
has been varied to obtain changing X values. The point 
at which all curves cross is at ^ = 1.0 and represents 
the afore-mentioned X^. Unlike the approximate Windenburg 

determination of this point, in this determination, Xj^ 



r 
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depends on the value of h/D as well as the fact that 
= 1. This dependence is borne out in the derivation 
of an expression for from formula #92A. This deri- 
vation is found in Appendix E and results in the 
equation : 



2.055Dp^ 

h 

which was used to plot h/Djjj versus 
values of Oy on figure 3. 




(A) 

E 

Xjjj for different 



Figure 3 may be used to determine for a given 
h/Dj^ and what X is necessary to prevent collapse be- 
tween the frames due to instability (\p <r 1.0). Areas 
to the right of the curve in question are areas of 
instability collapse, those to the left, primarily 
yield failure. As an example of this. Trilling's 
experimental results have again been plotted. Due to 
slightly varying geometries and Oy’s, all the models he 
used fall within the areas of the dots shown. The 
figure indicates that for Oy = 30,000 psi (see appendix) 
all the models failed at a below 1.0. Reference to 
figure 1 shows that this was indeed the case, and a 
further glance at Reference (19) indicates that a number 
of Trilling's models failed by lobar buckling or general 
instability, rather than by yield initiated collapse. 
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In Appendix A, a number of formulas and yield 
criterion are cited as being more accurate than the 
maximum circumferential stress criterion in determining 
collapse pressure due to yielding. One of these, Wenk's 
formula, has been investigated using the same models 
used in figures 1-3 in order to ascertain the influence 
of increased accuracy on the computation of Figure 

4 is the result of this investigation. 

Wenk applied formulas #92 and #92A to the Hencky- 
von Mises yield criterion at mid-thickness , mid-bay, 
and obtained collapse pressures within nine percent of 
actual experimental values. The Hencky-von Mises criterion 
in two dim.ensions is: 

Oy = - 0<pap + 

When it is applied to an unstiffened shell where 0^, = 

20p, yield failure occurs at 0^i = 1.160y, i.e., \p = 

1.16. Just as figure 2 indicated that for the maximum 
circumferential stress criterion, there is a X value 
where frame size has no effect on cylinder strength, figure 
4 shows that the same is true for the Hencky-von Mises 
criterion. However, the X value in this case occurs v;hen 
\p = 1.16 instead of 1.0. V^enk's formula: 
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1 



5F + .25 



like formula #92A, is actually the hoop stress relation 
times a factor (defined in Appendix E) , When it is solved 
for F by setting \p = 1.16, a Xj^ versus relation 

similar to formula (A) may be found. This formula is 
shown below and a detailed derivation of it is given in 
Appendix E: 



If plotted in the form of figure 3, this equation would 
result in a shift of all curves .IIX to the right of 
those obtained from #92A. This shift reflects the conser- 
vative nature of von Sanden and Gunther's formula 
and criterion. 

Thus far, equations and criteria have been developed 
to determine the extreme values of X between which 
stiffened cylindrical shell failure is primarily due to 
yielding of the material between the frames. Justification 
for using this as the design collapse mode comes from a 
number of sources. Wenk, a major advocate of the "one- 
hoss-shay" collapse philosophy, states that the approach 
of designing for yield failure is not inconsistent with 
that philosophy (221. Because of the effects of such 



X, 



m 




(B) 
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unknowns as material residual stresses and out of roundness, 
design for collapse by all three modes simultaneously 
is difficult, if not practically impossible. Therefore, 
the mode which is least affected by such strength modifiers 
should be the design collapse mode. Conversely, the 
collapse mode most affected by the modifiers should be the 
first mode taken into consideration and designed out of 
existence. Out of roundness is the primary reason for 
local failure of frames which can precipitate premature 
general instability collapse (15) . For this reason, 
collapse by general instability has been noted as being 
of prime concern in the preliminary design (10) . 

Reference to formulas for general instability given 
in the appendix will reveal that this collapse mode, unlike 
the others, is a function of cylinder length (between 
bulkheads), as well as frame area, moment of inertia, and 
spacing. Judicious selection of bulkhead spacing (if a 
free variable) and minimum frame sizes can eliminate 
this collapse mode from consideration. This should be 
done in the first stages of design development. In the 
final stages collapse pressure for this mode should again 
be checked and out of roundness taken into consideration. 

V7ith general instability designed for, collapse 
betv;een the frames due to buckling should be eliminated 




i 
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by insuring that the final design X will be between 0 
and Xjjj. This means that the yield strength of the material 
will be utilized as fully as possible by insuring against 
premature failure initiated by instability. In this way, 
the benefits of using high strength materials will not 
be nullified by instability collapse at low pressures. 

In the next section, considerations for the final 
step of designing for yield initiated failure between the 
frames will be developed. First, however, the question 
of frame adequacy must be looked into. 

Figures 1, 2, 3, and 4 were all developed from 
analysis of cylindrical shells with rectangular stiffeners. 
T-shaped stiffeners are in general use today however, so 
the question may be raised as to what effect the shape 
of the stiffeners has. An investigation of the formulas 
for yield failure previously cited will reveal that they 
are all based on the assumption of adequate frame strength. 
The only frame parameters employed in them are flange 
width in contact with the shell and' frame area. At X = 

Xj^, even these parameters cease to influence collapse 
pressure. Therefore, it can be stated that frame shape 
has no influence on formulas (A) or (B) . In general, 
as long as adequate strength is maintained, frame shape 
does not influence the collapse pressures of stiffened 
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cylindrical shells. 

The question of what constitutes adequate frame 
strength has been directly approached by a number of 
investigators (2,11,19). The standard practice of 
von Sanden and Gunther (20) and other early designers 
Vv'as to design the frame to be able to withstand 1.1 
times the total pressure on the length of plating it 
supported. This procedure, which employed the Foppl 
formula, was found to result in "excessively heavy frames 
(19) . In addition, the Foppl formula has been determined 
to give conservative results for frame buckling pressures 
In light of these facts, the frame design requirement 
chosen for the purposes of this report was that the frame 
be strong enough to support no more than the load on 
one effective length of shell. The conservativeness of 
the Foppl formula was accounted for by assuming that it 
predicted frames alternatively 4, 5, or 6 times too 
strong. The resulting formula used is; 

24EIj3 

P,., frame = p_ shell = X (4, 5, or 6) 

^ ^ 3t 

^b ^1 



Before proceeding into the actual design iteration 
resulting from the establishment of design criteria for 
stiffened cyl.ndrical pressure vessels under hydrostatic 
pressure, a recapitulation of the criteria is in order: 
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1. Upon definition of -operating requirements, initial 
Structure parameters should be established to 
minimize the possibility of general instability 
collapse of frame and shell. 

2. Possibility of collapse due to instability 
between frames should be reduced by limiting 
design choices to geometries located between 

m* 



X = 0 and X = X 
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III. VESSEL EFFICIENCY 



The measure of efficiency used in this paper 
was suggested by VIenk in (22) . It is slightly more 
complicated than the traditional Wei^ht/Displacement 
ratio in that it takes into account collapse pressure. 



10®W 



Pj, is the collapse pressure - Ibs/in^ 

V is the volume per unit length of cylinder - in^/ft length 

W is weight per unit length of cylinder - Ibs/ft length 



When used in a qualitative manner, this term clearly 
shows the relative efficiencies of different geometries. 
Hov/ever, it is not independent of material properties 
such as density and yield stress. The efficiencies in 
this report have been calculated using HY 80 as the 
material. In order to apply them to materials of other 
densities and yield strengths, they must be multiplied 
by the follov;ing factors; 






.284 

P 



For example, a geometry v/hich v/ould have an efficiency 

of 1.25 for HY 80, has an efficiency of: 

, . 50,000 ^.284 

1 . 25* r - - - - — * 



80,000 



.1 



.2.19 



- 22 - 



when the material is 50,000 psi aluminum. Once found in 
the above manner, the efficiency may be used to calculate 
the traditional Wt/Displacement ratio using the following 
conversion formula: 

Wt ^ 

A 3720n 

Reasons for this manipulation of n will become more clear 
at the end of the next section. 






1 
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IV. DESIGN ITERATION 



In the preceding section, design criteria have 
been developed which, if used, will produce a generalized 
cylinder geometry whose lvalue will be between 0 and 
This will result in a collapse pressure factor in excess 
of 1;0, thus insuring that collapse will be primarily 
due to yield and the yield strength of the material will 
be highly utilized. This section will outline development 
of an iterative method to find the particular X and ip 
values v.’hich result in maximum structural efficiency 
for a given collapse pressure. 

Figures 5, 6, and 7 are plots of efficiency against 
X for different values of h/Dj^^ and frame strength factors. 
As stated previously, Foppl's frame formula has been used 
in a modified form v;hich arbitrarily reduces its apparent 
conservativeness. In determining frame size for figure 

5, Foppl's formula was multiplied by a factor of 4 . In 
figure 6, the frame factor was 5 and in figure 7, it was 

6. The effect of this increasing theoretical strength 
of frames was to reduce frame size and increase apparent 
structure efficiency in succeeding figures. Since this 
is an arbitrary change, model tests are necessary to 
determine the actual frame factor which should be used. 

It was felt that a range of factors such as that covered 
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would indicate design trends which may be beneficial in 
general use without the necessity of establishing rigorous 
theoretical frame design criteria. Further substantiation 
of this method will be provided later in this section. 

The formula used to determine collapse pressures 
for these three figures was #92A (maximum circumferential' 
stress criterion) . Although this is now considered to 
be less accurate than a number of possible formulas, it 
has the advantages of always being conservative and 
having been widely used. The curves in figures 5, 6, and 
7 indicate that for a given h/Dj^j and frame factor, the 
efficiency increases with X until a maximum, is 

reached. Thereafter, the efficiency steadily decreases. 
v;ith increasing X . It is this point, hmax^ that the 
designer of stiffened cylindrical pressure vessels is 
interested in, therefore, only enough points have been 
plotted to indicate the hump of each efficiency curve. 

An interesting trend which is readily apparent in 
all three figures is that the efficiencies of all h/D 
combinations lie on the same line at the point where 
X = Xjj, for each h/D^^ (see figure 3). In figure 5, this 
line has a slope of about -.5. In subsequent figures, the 
slopes appear the same, so, since the lines begin at higher 
values, the final efficiencies will be higher. Of course. 
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extension of the lines upward to maximum efficiency also 
requires extension to the left to X = 0 and extension of 
the series of h/D curves to extremely thick cylinders. 

This extrapolation of computed data can be justified on 
the grounds that only very short, thick cylinders fail 
by pure yielding. Failure in this way takes the fullest 
possible advantage of material yield strength and results 
in the highest possible structural efficiency. Since 
efficiency in this report is also a function of displace- 
ment to weight ratio, the veracity of that previous 
statement, in the general sense, is dependent on material 
density. 

Another point which is illustrated in figures 5-7, 
is that for each h/Djjj there is only one value of X at 
which maximum efficincy is reached, regardless of frame 
factor. The effect of the frame factor is to change the 
actual maximum efficiency reached, but it has no effect 
on the X value where it is reached. This direct relation- 
ship between h/D^ and X at maximum efficiency is shown 
in figure 8. It may be explained in part by the fact that 
neither X nor h/Djj^ contain any frame size parameters and 
both Wenk's formula and formula #92A assume adequate 
frame strength. Also at that point of maximum efficiency, 
there are dis. rete values of ip which are dependent on 
frame factor. In order to show the relationship of frame 
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factor, A and at maximum efficiency, figures 9 and 
10 have been plotted. 

Before entering into the use of these figures, 
however, a slight digression is necessary. Recalling 
that formula #92A, used to obtain figures 5-7, was 
shown to yield less accurate collapse pressures than 
other formulas, it became necessary to investigate a 
formula of higher accuracy to determine if the same 
trends still held. The formula used was developed by Wenk 
and has been introduced in Section II. Data similar 
to that shovm in figures 5, 6, and 7, was developed for 
the same geometries using Wenk's formula. The actual 
results are given in Appendix G. Integrated curves from 
these results have been drawn on figures 8 and 10. 

Figure 8 is a graph of h/D versus A for points of 
maximum efficiency. The two curves represent results 
of the use of formula #92A and Wenk's formula. It will 
be noted that W^enk's form.ula shifts the curve to the right. 
Since this expression has been knov/n to predict collapse 
pressures in excess of actual values, it may be expected 
that the curve for actual structures will lie between the 
curves shown. No mention of frame size is made on this 
graph, as in neither formula did frame strength alter 
the relationship of h/D and A at maximum efficiency. 



nfi 
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Figures 9 and 10 are graphs of A. and iIj for 

#92A and Wenk's formula respectively. The three lines 
in each case indicate the three frame factors used. 
Comparison of these figures shows that increased accuracy 
does not change the trends previously noted. In fact, 
the n versus X curves are exactly alike. The 25 percent 
increase in efficiencies in figure 10 over figure 9 is 
another reflection of the conservativeness of formula 
#92A. 



With development of figures 8-10 complete, an 
iterative method to determine geometries of maximum 
efficiency using these figures can be outlined. Geometric 
parameters which can be found are h/Dn^ and X. Given 
material properties, distance between stiffeners may 
be computed. Actual stiffener parameters are left up to 
the discretion of the designer as the only requirement of 
the formulas used in this procedure is that the stiffeners 
be of adequate strength and stiffness. For comparison 
purposes, the frame parameters used in this report and 
those frame dimensional ratios coinciding with maximum 
efficiency structures are given in Appendix C. To begin 
the iterative procedure, the information v/hich must be 
specified is: 

1. type of analysis to be used 

a. maximum principle stress, #92A 
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b. Hencky-von Mises stress at mid-thickness, 
mid'-bay, #92 and #92A 

2. frame factor to be used: 4, 5, or 6 

3. ultimate collapse pressure desired 

4. material properties, Oy and E 
(for these curves E = 3 X 10^) 

5. outside diameter 

An example of the procedure is given on the next page. 
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Example 1: 
Given: 



Steps: 



Use of #92A, Frame Factor - 5 

Oy = 80,000 psi 

D = 40 ft. 

p^ = 1000 psi 



1. Initial h/Dj^ guess from hoop stress 



Pc 

(h/Dm) 

(h/Dn,) 

X 



2(h/D)Oy h/D = h/Dm 

1000/(2*80,000) 

.00692 

. 48 from figure 10 

1.29 from figure 11 
P 

c 




.*. P(, = 1290 psi - which is too high 



2 . 



Reduce h/Djn by division 



Pci 

(h/Dpfj) 

(h/Dm> 

X 



Pc 



1000 = 775 psi 

1.29 

775/(2*80,000) 

.00484 

.58 

1.41 

1.41*775 = 1090 psi - which is 
still too high 
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4 



Answer : 



Results ; 



. Reduce 



Pc2 

(h/Dn,) 

X 

Pc 



. Reduce 



Pc3 

(h/D^) 

X 

<P 

Pc 

X 

h/D 

n 

D 

h 

Wt. 



h/Djj^ again 
775 

= T7o9- " 

= 710/(2*80,000) = .00443 

= .606 
= 1.44 

= 1.44*710 = 1020 psi - which 

is close enough, but try a 4th 
iteration 



h/Dm 



710 

1.02 



697 psi 



= 697/(2*80,000) = .00436 

= .613 

= 1.45 

= 1.45*697 = 1012 psi 



= .613 

= .00436 

= 1.22 

= 1.45 

= 40 ft. 

= 2.1 in. 

• ' 1.56 “ft. 

= 18,000 Ib/ft 



conventional V7t/Disp 



.222 
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V. DISCUSSION OF RESULTS 



The results of this report may be divided into two 
classes; those due to analysis and those due to synthesis. 
Analysis by formula #92A and Wenk's Formula at low 
lambda values has shown that considerable gains in 
efficiency can be realized by placing frames close together 
and reducing shell thickness. This results in 4 > values 
greater than 1 and helps to realize the stress potentials 
of high strength materials. The extrapolation of the 
design curves in figures 11 and 12 appears to justify 
the use of sandwich construction for very high pressures. 
Such construction is analogous to design in the very low 
lambda region, and a discussion of this is included in the 
appendix. Inferences which may be drawn from the treat- 
ment of very lov; lambda design in Appendix B concern the 
apparent inapplicability of X and in this region. This 
results from the fact that in this region, frame strength 
contributes so much to overall strength that it cannot 
be ignored as it is in both X and ip . 

At the high end of the low lambda region, maximum 
structural efficiency can be realized only in relatively 
low pressure vessels. This is shown in two examples in 
the appendix. 
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The result of synthesis in this report is the 
iterative method outlined in Section IV. By its very 
nature, this method results in thin-skinned vessels 
with frames relatively close together. The questions 
which this geometry may raise concerning modes of collapse 
other than yield-initiated, must be answered by indepen- 
dent analyses. It will be noted, however, that the 
geometries used to plot the cuifves of figures 8-13 were 
analyzed for failure both by instability between frames 
and general instability and the indicated collapse mode 
for all was yield-initiated failure between rings. 
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VI. CONCLUSIONS 

1. The dimensionless parameter, X, may be used to 
define all stiffened cylindrical pressure vessel geometries 
except those in which the frames provide little support 
because of their wide spacing, and those vessels in which 
the frames provide most of the strength by virtue of their 
size and close spacing. 

2. 7^bsolute structure efficiencies in the low lambda 
region are highly dependent on frame spacing and only 
slightly dependent on frame size (as long as the frames 
are of adequate strength) . 

3. Relative structure efficiencies in the low lambda 
region are virtually independent of frame size. The lambda 
values at which maximum efficiency is reached depend on 

the h/Dj^ ratio and are not affected by changes in frame size. 

4. Structures with values well in excess of 1.0 
are reasonable and generally provide increased efficiencies 
over those where is limited to 1.0. 

5. As the design collapse pressure increases, the 
value of lambda which coincides with maximum structural 
efficiency decreases. This implies that higher efficiencies 
can be obtained by reducing shell thickness v/ell below that 
required by the hoop stress relation and simultaneously 
increasing the number of frames. 

6. At extremely low lambda values, general instability 
becomes the most likely mode of failure. 
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VII. RECOMI-IENDATIONS 

Since this thesis was limited to studying the 
parameter, X, and its effect on efficiency, other 
definitive parameters, such as those nov; used at the 
David Taylor Model Basin should be investigated for 
their effect on efficiency. It has been noted that the 
lambda term is relatively useless when the frames are 
very close together and are very strong. This condition 
should be corrected by defining a new parameter similar 
to lambda in its application, but allowing for the effect 
of frame strength in geometric comparisons. 

The basic approach of this thesis should be used 
to extend the curves of maximum efficiencies, found in 
Section IV, into the regions of higher collpase pressures. 
With this extention, a thorough re-examination of general 
instability collapse in these regions should be undertaken. 
Finally, experimental tests of maximum efficiency geo- 
metries, as developed by the design iteration, should be 
undertaken to provide evidence of the method's worth. 
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APPENDIX A 

HISTORICAL BACKGROUND 

Prior to 1917, the design of stiffened cylindrical 
pressure vessels was based on the assumption that vessel 
strength depended primarily on frame strength. The shell 
was treated as a flange of the frame and strength was 
calculated from equations which defined collapse by 
instability of frame rings under radial loads. A widely 
used formula for this was developed by Levy, but is 
commonly called Foppl's Formula (19,24): 




While this formula was relatively accurate for the 
time, and where frames alone were involved, when it was 
applied to a stiffened cylindrical shell, its shortcomings 
became apparent as higher collapse pressures were required. 

In order to improve upon this frame-oriented theory, 
von Sanden and Gunther derived expressions based on the 
observation that stiffened cylindrical shells were not 
collections of connected rings, but, in fact, were closed 
cylinders reinforced by ring frames. Deciding that the 
shell was the real strength member, they predicted failure 
due to shell yielding at the area of maximum stress. This 
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ma:timum principle stress criterion could be satisfied 
by the longitudinal stress at the frames or by the 
circumferential stress midv;ay between them. Formulas 
for both stresses were developed and the well-known 
formulas #92 and #92A resulted when pressures were 
for in terms of the stresses. Substituting material 
yield stresses into #92 and #92A resulted in two yield 
pressures for the structure, of which the lower one was 
said to predict failure. It must be emphasized that 
failure was predicted on the basis of the onset of 
yielding. 

For a number of years after the 1918 publication of 
these formulas, #92 was thought to give the critical 
collapse pressure since the greater stress was longitudinal. 
In 1935, however, this assumption was found to be in 
error when circumferential stress was determined to be 
critical in the collapse of thin-shelled pressure vessels 
(19) . 

The von Sanden and Gunther formulas are: 



#92 





h+ 1.8 IK 



85 - B 

1 + B , 
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#92A p. 



y 




h + H (.85 - B\ 

\ 1 + B/ 



The variables K, H, B, and 6 are defined in the Formula 
section. 

About this same time, it was generally recognized 
that the area between stiffeners very rarely collapsed 
as a result of pure yielding. In anything but thick 
cylinders, collapse occurred from a combination of 
yielding and buckling v;hen the distance between stiffeners 
(effective length) was below a critical length. This 
critical length was dependent on shell L/D and h/D ratios 
as well as material properties. When the distance between 
stiffeners was above the critical length, the tube was 
considered to be of infinite length and the formula of 
Breese and Bryan for collapse by instability could be 
applied: 



2E 



1 - 




For cylinders below the critical length, formulas 
for collapse by instability were developed by von Mises 
and Southwell. Windenburg simplified von Mises lengthy 
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forraula and, by sacrificing . a very small degree of accuracy, 
came up with the following formula (25) : 



2.42E (h/D) 

(1 - y V [L/Dro - .45(h/D^) ] ‘ 



At the time Windenburg developed the above formula, 
a Pc versus L/D coordinate system was commonly used to 
plot experimental and theoretical results for changing 
values of h/D. In an effort to develop an analog to the 
slenderness ratio, L/r, in column theory, and aid comparison 
of results, he developed a nev7 coordinate system based 
on the coordinates: 






2h/Dm 






(L/Di!i) 



(lOOh/Dro) 



Windenburg then corrected his coordinates for the physical 
properties of materials in the following manner and came 
up with the well-knov/n variabels, slenderness ratio X , 
and pressure factor, ip . 



2h/Dj(jay 
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\ 



1000a, 




L 




- <VDm) ^ 

(h/D„) > 




This done, Windenburg rewrote his instability formula in 
terms of his new variables: 



'P 



1.30 



X ^-e 



where e 



.045 



lOOOOy / 1 




1 



E \ .91 / lOOt/D 



may be neglected in general to get 






1.30 






This formula is plotted on figure A-1 and may be 
compared there with the plot of von Mises' formula and 
numerous experimental results. One thing which may be 
seen at this time is the good agreement these formulas 
have with eachother and the experimental results above 
a X value of 1.4. 

After a decade of relative inactivity, the investi- 
gation of stiffened cylindrical pressure vessels was renewed 
in the early 1950 's. The main intent of the investigations 
of this period was to correct the number of theoretical 
deficiencies in formulas #92 and #92A. The first two 
deficiencies were that the formulas predicted linear 



1 ’ 
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relationships betv/een pressures and stresses and ignored 
the effect of moments developed due to pressure at the 
ends of closed cylinders. Both problems were corrected 
by Salerno and Pulos (17) , who developed stress formulas 
similar to #92 and #92A but taking both the above 
moments and the influence of axial stress in the shell 
tending to expand the frame into account. Their analysis 
resulted in a non-linear relationship between pressures 
and stresses which, however, did not alter the results 
of stress calculations significantly from those found by 
#92 and #92A. 

The next problem found with von Sanden and Gunther's 
theories had to do with the inadequacy of the maximum 
principle stress criterion. Use of the Hencky-von Mises- 
yield criterion: 

^ - OtOp + 

at the exterior mid-bay point resulted in predictions closer 
to actual collapse pressure than obtained with the former 
criterion . (9) . However, here again collapse pressures 
were associated with the onset of yielding and no allowance 
for plastic reserve strength was made. 

This condition was corrected by Wenk et. al. (21) in 
a somewhat haphazard manner by taking the Hencky-von Mises 
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yield criterion at mid-bay and mid-thickness . Although 
Wenk termed this a "necessary, if not sufficient empirical 
yield condition to produce failure", experimental results 
have since indicated that thought this formula predicts 
within nine percent of actual collapse pressures, it 
tends to predict high. (9) . An interesting observation 
on the use of the Hencky-von Mises criterion is that for 
an unstiffened cylindrical shel'l, Ot = 20p. This results 
in failure of such a shell occurring at a«p = 1.160y; 
stated in Windenburg's variables, = 1.16. 

The use of form.ulas #92 and #92A in the development 
of Wenk's formula appears to be the last time these 
long-standing elastic stress determinators were used in 
a yield criterion. The elastic stress formulas of 
Salerno and Pulos (17) , having been recognized as more 
rigorous and slightly more accurate, were used in the next 
collapse criterion developed by Lunchick. Lunchick 
utilized the Hencky-von Mises yield criterion, however, 
he specified that collapse did not occur until the load- 
carrying capacity of the shell v/as completely exhausted 
at the critical mid-bay point. Lunchick 's three-hinge 
mechanism stemmed from his assumption of plastic hinges 
forming at the frame at the same time a completely plastic 
hinge had developed at mid-bay. Lunchick 's method predicted 
collapse pressures with better accuracy than any other 
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in the range of low X values. However, its inclusion of 
design curves to determine cylinder plastic reserve 
strength has made it more difficult to apply in computer- 
based design methods. 

It is now generally accepted that stiffened cylindrical 
shells fail by one of three distinct modes or any 
combination of them. These modes are (14) : 

Mode 1. Axisymmetric collapse of the shell between 
adjacent ring frames. 

Mode 2. Asymmetric collapse of the shell between 
adjacent ring frames (lobar) . 

Mode 3. Overall asymmetric collapse of the shell 
and frames together. 

Mode one had traditionally been referred to as 
axisymmetric shell yielding . (9 , 22 , 23) . It has been 
pointed out, however, that this mode is the result of a 
combination of yielding and axisymmetric buckling (14,22) 
and the main difference between this mode and mode two is 
the difference in resulting collapse patterns, not basic 
collapse mechanisms. Axisymmetric collapse pressures 
have been calculated with good accuracy through use of 
yield criterion alone though, and it is recognized that 
pre-collapse stresses, as determined by von Sanden and 
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Gunter (20) or Salerno and Pulos (17), are the primary 
cause of failure in this mode. 

The problem of axisymmetric collapse by elastic 
buclcling has been solved by Salerno and Pulos in (13) . 
Lunchiclc (10) has investigated axisymmetric failure by 
inelastic buclcling and derived the following expression 
for it: 



von Mises' and Windenburg ' s formulas for the calcula- 
tion of failure pressure by elastic asymmetric (lobar) 
buckling betv;een ring frames have been given. A more recent 
solution to this problem has been found by Reynolds and 
is outlined in reference (14) . Reynolds has also developed 
a formula for the inelastic buckling of stiffened cylindri- 
cal shells which appears to be reasonably accurate, but 
is consistently high in its calculation of collapse pressures 



Per 





ir^Rh 



( 1 ) 



Mode thr e, collapse by general instability, has 
been found to be of primary improtance in the preliminary 
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design of pressure vessels (14). Solution of elastic 
general instability collapse was first undertaken by 
Tokugawa, who was also the first to identify this mode. 
His formula for collapse pressure, (which was later 
confirmed by Bryant) , is as follows: 



Per 



Eh 






_n^ - 1 + j ) (n^ + dM ^ 



(n^ - 1) 



Ele 






where X = — 

^b 

The original Tokugawa formula as given in (19) has a 
number of extra terms in it which are seen to be 
negligible for geometries of general interest. 



The lengthy formulas developed by Kendrick for 
elastic general instability are given in (14), however, 
a graphical solution of Kendrick's equation as developed 
by Reynolds, and refined by Ball, is included as figures 
A-2 and A-3. In this solution, the sum of two effects, 
each represented by a graph, is used to determine collapse 
pressure. Figure 2 is the shell effect and figure 3 
the frame effect which reflects the bending stiffness of 
the cylinder. Figures 2 and 3 also indicate the Bryant- 
Tokugav/a solution. Collapse pressure by either set of 
curves is found in the following manner: 
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FIGURE A -2 



(SHELL EFFECT) 
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FIGURE A -3 (FRAl'lE EFFECT) 
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In order to round out this development of stiffened 

cylindrical pressure vessel collapse theories, mention 

must be made of the effect of actual physical manufacturing 

conditions. In all but the sma'llest and most painstakingly 

fabricated cylinders, imperfect production processes are 

the cause of a number of collapse strength modifiers. 

Those modifiers introduced into the material itself are 

in the form of stresses due to fabrication, rolling, and 

v/elding of flat plating. In thick homogeneous plating, 

the homogeneity itself may be questionable and thus cause 

stress concentrations. A major strength modifier is the 

imperfect circularity of most welded cylinders. This 

modifier has received some theoretical consideration as 

✓ 

it has a considerable effect on the resistance of a stiff- 
ened cylinder to failure by general instability. A 
formula for this is due in part to Kendrick (2) : 




.85 




'^max 



1+6 



As long as the a 



max 



is below material yield stress 



collapse by general instability is not indicated. 
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This section has been included to give a brief 
background of the development of theoretical collapse 
predictors for stiffened cylindrical shells. In particu- 
lar, a history of the variables of X and ip is helpful 
in understanding their use in this paper. 
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APPENDIX B 

DISCUSSION OF VERY LOW LAMBDA DESIGN 

As noted in Section IV, overall structural efficiency 
approaches its highest value as rp goes to 1.0 and X goes 
to 0. For T-shaped frames of the geometry used in this 
report, X values in the area below .2 may be considered 
very low lambda design. For other frame geometries, this 
may be higher. For example, the sandwich construction 
in figure A-4 has a X value of about .337. The frames used 
there have thicker webs than flanges, whereas the reverse 
is true for the frames of this report. The membrane shell 
design, and other complications resulting from the use of 
high strength titanium also plays a part in deciding the 
optimum X value. 

It was noted in Section IV that the danger of general 
instability collapse would be greatly increased by the 
use of a very low lambda design. This has been shown to 
be true in model tests of cylindrical shells with 
geometries similar to that of figure A-4. In fact, in the 
testing of these models, a new mode of collapse was observed 
- that of plastic general instability - and new theories 
had to be developed to predict it. It may also be observed 
that collapse of these models took place at \p value in 
excess of 3.0. This suggests that at very low lambda 
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values, , becomes a rather meaningless term unless 
the framing is very small in comparison to the shell 
thickness. Indeed, the large effect of framing at these 
X values throws into doubt the meaningfulness of X itself. 
This doubt has been substantiated but not fully explained 
in reference (6) . 
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APPENDIX C 
FRAME PARAMETERS 



A T-shaped frame was used in the calculations 
leading to figures 5-10. The proportions of this 
frame were taken from a work by McGinley (11) and are 
as follows: 

frame width/frame depth = 1.0 

frame width/web thickness = 8.5 

flange thickness/web thickness = 1.7 

Figure A- 5 is a graph of b/D versus X for the frames 
for both formulas #92A and Wenk.'s., Since the geometry 
of optimized vessels was found to be a function of required 
collapse pressure, and not required diameter, the optimized 
frame dimensions should be proportional to the major 
cylinder dimensions., hence the b/D parameter. Figure A-4 
may be used in conjunction with figures 8, 9, and 10 to 
determine the frame dimensions. Care should be taken, 
however, to match yield formulas and frame factors when 
using these figures together. 



- 64 - 



APPENDIX D 



EXAJ:4PLES of design iteration 
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Example 2: 

This example is included to indicate the results 
of the iterative method when applied to relatively low 
pressure vessels. 



Given: Use of #92A Frame Factor 5 

Oy = 30,000 psi 

P(^ = 400 psi 

D = 40 ft. 



Steps : 



1. Initial h/D^^j guess from hoop stress 



(h/D^) 

(h/D^) 

X 

'I' 



Pc 



400/(2*30,000) 

.00666 

.495 

1.30 

1.30*400 = 520 psi 

which is too high 



2. Reduce h/D^p by division 



Pci 



(h/D^) 

X 

'I' 



400 



= 308 psi 



1.3 

= .00512 

= .565 

= 1.39 

= 1.39*308 = 428 psi 

v/hich is still too 






act 



428 



high 
= 1.07 



400 
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3 . Reduce h/Dj^ again 



Answer ; 



Results: 



Pc 2 

(h/D^) 

X 

'P 

Pc 



h/Dn, 

m 

ti 

n 



308 = 288 psi 

1.07 

.00480 

.584 

1.41 

1.41*288 = 406 psi 

which is close enough 

.584 

.00480 

1.27*30,000 

80,000 

.476 

1.41 



D 

h 

L 

Wt 

A 



= 40 ft. 

= 2.3 in. 

= 4.5 ft. 

= ^ 

3720n 



406 

3720(.476) 



Wt 

A 



.23 
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Example 3: 

This example is exactly like example 2 except that 
a higher strength material is used. 

Given; Use of #92A Frame Factor 5 



Q 

11 


50,000 psi 


Pc = 


400 psi 


D 


40 ft. 


Answer; X 


= .825 


h/Dxn 


= .00232 


n 


= 1.11* 50,000 

80,000 


n 


= .694 


'P 


= 1.735 


Results; D 


= 40 ft. 


h 


= 1.11 in. 


L. 


= 1.83 ft. 


Wt 

A 


= .157 




I 

I 



1 
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Exaraple 4: 

This example is included to indicate the use of the 
frame parameter curves in Appendix C. 

Given: Use of #92A Frame Factor 6 



= 


100,000 psi 


D 


40 ft. 


Pc = 


1500 psi 


Steps: 





1. Initial h/D^j^ guess 



(VDn,) 


= 1500/(2*100,000) 


(h/Dn,) 


= .007 


X 


= .48 


'P 


= 1.33 


Pc 


= 1.33*1500 = 2000 psi 



2 . Reduce h/Dj^ 



Pci 


= 1500 _ 

1.33 


(h/Dj[j) 


= 1130/(2*100,000) 


(n/Dm) 


= .0056 


X 


= .54 


'P 


= 1.40 


Pc 


= 1.40*1130 = 1580 psi 


1'aot 


- 1580 _ 1 05 

1500 
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3. 



Reduce h/Dj^- again 



Pc2 = 


1130 

1.05 


= 1075 psi 


(h/D^) = 


1075/(2*100 


,000) 


(h/D^) = 


.00537 




X 


.55 




iP 


1.41 




Pc 


1.41*1075 


= 1515 psi 




- 


close enough 



Answer: 



Results : 



X = .55 

(h/Dj^) = .00537 



100,000 

80,000 

ip = 1.41 

^ = .00248 

D =40 ft. 

h =2.6 £n. 

L = 1.42 ft. 



Vlt 

A 



.252 



Frame Dimensions: 
b 
w 
f 

t'2 



1.18 in. 

8.05 in. 
10.1 in. 

2.05 in. 



1.59 
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APPENDIX E 

FORMULAS AND DERIVATIONS 



In this appendix, the formulas used in the computer 
programs in Appendix F are listed. Also, the derivation 
of formulas (A) and (B) , Section II, and the solution 
of the Hencky-von Mises criterion for unstiffened tubes 
are presented. 



FRAI4E FORMULAS: The formulas below are used to 

predict frame failure by instability. Collapse of the 
frame by yielding v;as not investigated as the frames 
were not required to hold up after shell failure and the 
maximum stress in the frame is always less than the 
maximum stress in the shell (19) . However, this stress 
relationship may not hold in the very low lambda region. 



Foppl Formula: 
Pc = 



24EIj3 

Db'Li 



Formula #88 

von Sanden and Gunther: 



Pc 



24E1]^ (1+B) 
Db^b (1 + 



xs 



bh 



where B 



A + bh 
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0 



1.82 



L/Dm 

h/Dn, 



6 



L NB 
b 0 



N = cosh0 - sin0 
sinh0 + sin0 



SHELL FORMULAS: The first formula below predicts 
collapse of the shell between frames due to instability. 
The next four predict collapse due to yielding of the 
material according to different yield criteria. 

Formula #92 

von Sanden and Gunther 

Failure due to axial stress exceeding yield stress: 
n ha / 1 \ 

Dm y .5 + 1.81KC 1 

_ sinh0 - sin0 
sinh0 + sin0 

o _ .85 - B 
1 + 6 

B and 6 are defined in Formula #88 



Formula #92A 

von Sanden and Gunther 

Failure due to circumferential stress exceeding yield 
stress : 



Pc = 



( ^ ) 

' 1 + HC / 
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H = -2 



(1 + M) sinh^ (1-M) coshy 



sinh0 + sin0 



M = 



3v' 



1 - u" 

C is defined in Formula #92 



Wenk ' s Formula 

This formula was actually developed by Wenk and others 
at DTMB using #92 and #92A to determine the principle 
stresses. These were applied at mid-bay and mid-thickness 
using the Hencky-von Mises yield criterion. 






F^ - .5F + .25 



F = 1 + 2QC 



Q = - 



. , 0 0 ^ , 0.0 

sinh'j cos-^ + coshj sin^ 



sinh0 + sin0 



0 and C are defined in Formula #92 



Lunchick's Formula 

This formula uses the principle stress relations developed 
by Salerno and Pulos rather than #92 and #92A (9) . It 
applies them to the mid-bay location using the Hencky- 
von Mises criterion, however, failure is not predicted 
until the yielding which occurs first in the outer surface 
develops into a fully plastic hinge. Lunchick's formula 
predicts this yield pressure, but curves (included here 
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as figure A-6) must be used 'to determine the actual 
collapse pressure. 



P 



c 



oJLt 

m 




_JL 

XZ + 



X = 1 + HC 
Z = .5 + 3.62TC 

. , e e , e . e 

sinhj cos-^ - ,cosh^ 
sinh6 + sin6 

6, C, and H are defined in Formula #92’A 



Entry into the curves requires computation of the 
following ratios: 

Circumferential Bending Stress _ ^BT 
Circumferential Membrane Stress 

Longitudinal Membrane Stress _ 

Circumferential Membrane Stress 



^m Pc 

Op - 

4h 

0^ = 2ap(l + 2QC) 

— 2 Op TC 



Collapse Pressure 
Yield Pressure 
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Curves for Lunchick's Collapse Pressure 
Factor - From DTMB Report 1291 



FIGURE A -6 



Longitudinal Membrane Stress 
Circumferential Membrane Stress 
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GENERAL INSTABILITY FORMULAS: While Tokugawa ' s 

formula and graphs of Kendrick's third solution for 
general instability have been given in Appendix A, 
Bryant's simplification of Kendrick's first solution 
has been used in the computer program. 



Bryant's Formula: 



cr 



2 ^ E 
m 



(n^ - 1 + ) (n^ + X^)' 



+ 



n 



2 



1 



R 



m 



-Ell 



X = an integer 

^b 

n = the integer which gives the lowest Pcr* 



Since out of roundness effects general instability 
collapse, a formula taken from (2) has been used in the 
following form: 



'^max 




1 -I- 



1 -I- 



.853 

_B 

3 
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e = eccentricity 

p^= desired collapse pressure 

n, Per' ® determined in previous 

formulas . 
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FROM PRINCIPLE STRESS CRITERION AT MID-BAY EXTERIOR, 
in ~ — — 

FORMULA #92A 



1. Formula #92A 



o hqy 

_Dhl_ 



1 + H 



.85 - B 
1 + 3 



2. ij; 



> h 0 

•BiT y 






1 @ 



1 + H 



.85 - B 
1 + 6 



* 1 



1 + H 



.85 - B 
"1 + 6 



= 1 



H 



\85 - B 

V 1 + 8 



= 0 



H = 0 



H = - 



A 



e 



3 sinhj cosj + cosh^ 



sinh + sin 

H=O@0 =3.74 radians 

0 = (1.82 L/Drn)/\|l^ 

L/Dnj= 2.055 \[1^ 



■"m 



m 



h/D^= .236 (L/Dj^) 
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3. X = 



4 


1 


Oy 


\ 


<h/Dm) ’ 1 


E 




2,055 Dyn 


0y 


\ 


h 


E 



Equation (A) , Section II 
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fro m HENCKY-VON MISE S CR ITERION AT MID-BAY AND MID- 
T HIC KNESS , 

VON SANDEN AND GUNT H ER FORMU LAS 



1. Wenk's Formula 
2 hcFv 



-'m 



- F + 1 



F 


= 


1 + 2Q 


/ . 


85 - B \ 








1 


1 + B / 








9 


i i 


Q 


= 


sinh 


2 


cos 2 + cosh 2 sin 



sinh 6 + sin 6 



2. \p 






9 



A 



1.1540 ^ 



F^ - — + i 

2 4 



F^ - 



H + i 

2 4 



•= 1.154 



.75 



F = .25+ .75 



.*. 1 = 1 + 2Q / .85 - B \ 

V 1 + e / 



. Q =009=4. 68 radians 



(1.82 L/D^) / \[h/D^ 
L/D^ = 2.565 fh/Dj^ 

h/Dj^ = 1 152 (L/Dj^) 



2 
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3. X = “ 


(L/D^) ^ 

_ 1 




\ 


(h/D^) ^ \ 


E 




2.565 Dj^ 




h 


E • 



Equation (B) , Section II 
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SHELL ? ^ JIDELY SPACED STIFFENER S 

1. He' --von Mises Criterion 

'^3 ~ three principle stresses 

2c,. = (C]^ -02 + (^2 -Oj )^ + (03 - 



For 



unstiffened tube; 



2a- 



2c^/ = 02^+02^ + 402^ 

Oy“ -■= 3o| 



0 



o 

2 



o 

1 



1 

1.732 Oy 
2 



1.732 



.5770y 

1.1540y 




I 
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APPENDIX F 
COMPUTER, PROGRAMS 



Two computer programs were used to produce theoretical 
data for this report. Program 1 utilizes equation (A) 
of Section II to determine curves of h/Dj^ versus 
for eight values of Cy. Since equation (A) is based on 
the 0 value at which the variable H in formula #92A 
goes to zero, the program first solves for that 6 and 
then goes on to vary h/Dj^ and Oy, solving for This 

program was used only to define the upper limits of the 
yield failure area and so, is not actually necessary for 
the iterative optimization proceedure of Section IV. A 
flow chart and variable chart are included with the program 
listing and printout. 

Program 2 follows a scatter-gun philosophy in that 
it determines collapse pressures of given geometries for 
all three modes of failure by various formulas. All of 
the formulas are mentioned in Appendix A and completely 
defined in Appendix E. The input of this program consists 
of outside diameter, desired range of collapse pressures, 
frame proportions, material properties, and range of 
values. Shell thickness is determined by the hoop stress 
relation and frame size by a modified Foppl formula. Program 
2 was not intended to be used to produce a complete 
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optimized design for a given set of requirements. Its 
main use, in solving for and printing collapse pressures 
for all modes, is pointing up trends and areas where 
various formulas fail. It has provisions for any 
number of solutions to be printed and any geometry to 
be investigated. In short, it is meant to be used as 
a tool by a designer who optimizes, using judgment and 
experience as well as computer time. A flow chart and 
variable chart are included with the program listing 
and a sample of its printout. 



o o o o 
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PROGRAM 1 



THIS PROGRAM DETERMINES LAMBDA MAX FOR CHANGING H/D 
AND YIELD STRESS USING A FORMULA DERIVED FROM 92a . 



P2 = .09 
NUM = 1 
NEND = 50 
TH = 3*7 
DTHP = .1 
DTHM = .1 
1 NUM = NUM + 1 

IF(NUM.EO.NEND) GO TO 30 
THT = TH/2. 

DEN = SINH(TH)+SIN(TH) 

EHO = (l.+(3*P2/(l-P2)) ^ 5 INH (THT)*C0S (THT) 

EHl = (1 .-(3*P2/(1-P2))**.5)*C0SH(THT)*SIN(THT) 

EH = -2-^(EH0+EH1)/DEN 
IF(EH) 19,30,9 
9 NP = NP+1 

IF(NP.EQ.3) GO TO 10 
DTHP = DTHP/2 
10 TH = TH-DTHM 
M = 1 
GO TO 1 

19 M = M+1 
1F(M.EQ.3) GO TO 20 
DTHM = DTHM/2 

20 TH = TH + DTHP 
NP = 1 

GO TO 1 

30 WRITE(6,180) NUM,EH,TH 

180 FORIATC AFTER’, 13,' ITERATIONS, H IS',F8.6,' AND THETA 
E = 30000000.0 IS' F6 3) 

DOH =0.0 * • / 

DDOH =40.0 
NDOH = 15 
DY = 10000.0 
NY = 8 

DO 50 N = l.NDOH 
WRITE (6, 2 05) DOH 
Y = 30000.0 
DO 40 M = 1 ,NY 
ELAM = (TH/1 .82*DOH*Y/E)**.5 
WRITE(6,210) Y,ELAM 
40 Y = Y+DY 
50 DOH = DOH.+ DDOH 

205 FORIATC U. ©DAS FOR DM/H = ',F7.2,' ARE ') 

210 F0PJAT(2F12.4) 

END 




I 
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FR INTOUT OF FRUG f L^K 1 (Sample) 



AFTER 50 ITERATIONS, H 
LAMBDAS FOR DM/H = 



IS -.000000 AND THETA IS 3.74 
0.0 ARE 



30000.0000 


0.0 


40000.0000 


0.0 


50000.0000 


0.0 


60000.0000 


0.0 


70000.0000 


0.0 


80000.0000 


0.0 


90000.0000 


0.0 


100000.0000 


0.0 


LAMBDAS FOR DM/H 


= 40. 


30000.0000 


0.2867 


40000.0000 


0.3311 


50000.0000 


0.3702 


60000.0000 


0.4055 


70000.0000 


0.4380 


80000,0000 


0.4683 


90000.0000 


0.4967 


100000.0000 


0.5235 


LAMBDAS FOR DM/H 


= 80.1 


30000.0000 


0.4055 


40000.0000 


0.4683 


50000.0000 


0.5235 


60000.0000 


0.5735 


70000.0000 


0.6194 


80000.0000 


0.6622 


90000.0000 


0.7024 


100000.0000 


0.7404 


LAMBDAS FOR DM/H 


= 120.1 


30000.0000 


0.4967 


40000.0000 


0.5735 


50000.0000 


0.6412 


60000.0000 


0.7024 


70000.0000 


0.7587 


80000.0000 


0.8110 


90000.0000 


0.8602 


100000.0000 


0.9068 


LAMBDAS FOR DM/H 


= 160. ( 


30000.0000 


0.5735 


40000.0000 


0.6622 


50000.0000 


0.7404 


60000.0000 


0.8110 


70000.0000 


0.8760 


80000.0000 


0.9365 


90000.0000 


0.9933 


100000.0000 


1 .0470 
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f START 
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NOMEMCLATURE OF PROGRAM 1 



DOH 

DDOH 

NDOH 

E 

EH 

ELAM 

P2 

TH 

Y 

DY 

NY 



Incremental change in Dm/h 
Number of changes in Dm/h 
Young ' s Modulus 

Variable H of yield failure formulas 
^m 

Poisson's Ratio squared 

Variable 9 of yield failure formulas 

Material Uy 

Incremental change .in Uy 
Number of changes in Oy 



PROGRAM 2 



THIS PROGRAM DETERMIR’ES THEORETICAL COLLAPSE PRESSURES 
FOR STIFFENED CYLINDRICAL PRESSURE VESSELS . THREE MODES 
OF COLLAPSE ARE INVESTIGATED WITH EMPHASIS ON YIELD 
PRECIPITATED BUCKLING BETWEEN THE FRAMES. FRAME SIZE 
IS DETERMINED USING A MODIFIED FOPPL FORMULA. ALL INPUTS 
ARE HANDLED BY MAIN PROGRAM CARDS - NO DATA CARDS ARE 
USED. 

MF = 6 
D = 480.0 
PC = 300.0 
DO 700 NU = 1,14 
Y = 80000.0 
E = 30000000.0 
R1 = 1.0 
R2 = .11764 
R3 = 1.7 
P = .3 
P 2 = P ^^*2 
SR = .38 

1F(PC.LT. 1000.0) SR s .46 
DSR = .02 
NSR = 20 
T = D*PC/(2*Y) 

DM = D-T 

TOD = T/DM 

WRITE (6, 160) T,TOD 

DO 500 IN = 1,NSR 

EL = SR**2*DM*E/Y*(T/DM)**1 .5 

FRAMDl = 12.0 

DFRP =1.0 

DFPJ-I =1.0 

NUM = 1 

NEND = 70 

RATIO = £L*FC/(24*E*mf) 

1 NUM = NUM + 1 

IF(NUM.SQ.NEND) GO TO 30 

F = FPu\MDl*Rl 

B = FRAMDl *R2 

T2 = B*R3 

W = FRAMDl - T2 

YB = YENTRD(B,T,B,W,T2,F) 

YBi = yi:oh(b,t,b.w,t2,f,yb) 

DBT = (D-2*(YB+T))**3 
RAT = YBI /DBT 
IF(RAT-RATD) 19,30,9 
9 NP = NP+1 
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IF(NP.EQ.3) GO TO 10 
DFRF = DFRF/2 
10 FRAIIDl = FRAMDl - DFRM 
M = 1 
GO TO 1 

19 M = M+1 
IF(M.EQ.3) GO TO 20 
DFRM = DFRM/2 

20 FRAMDl = FRAMDl + DFRP 
NP = 1 

GO TO 1 

30 YF = YENTRD(0.0,0.0,B,W,T2,F) 

Y1 = YENTRD(EL,T,B,W,T2,F) 

YI = YMOM(0.0,0.0,B,W,T2,F,YF) 

Yll = YII0M(EL,T,B,W,T2,F,Y1) 

YBIOYI = YBI/YI 
YIIOYI = YII/YI 
RONE = (D-2*(T+Y1))/2 
DB = D-2*(T+YB) 

A = F*T2 + B’'‘-W 
AB = A+B*T 
EB = B*T/AB 

THETA = 10*(12-12*P2)**.25 
TH = THETA*EL/DM/(100*T/DM)**.5 
THT = TH/2 

DEN = SINH(TH)+SIN(TH) 

EN = (COSH (TH) -COS (TH)) /DEN 
EK = (SINH(TH)-SIN(TH))/DEN 
BETA = 2 .*EL*EN*EB/(B*TH) 

EHO = (1+(3*F2/(1-P2))**.5)*SINH(THT)*C0S(THT) 

EHl = (1-(3*P2/(1-P2))**.5 )*coSH(THT)*SIN(THT) 

EH = -2* ( EH 0+ EHl ) /DEN 

EQ = -(SINH(THT)*COS(THT)+COSH(THT)*SIN(THT))/DEN 
ET = -(SINH(THT)*COS(THT)-COSH(THT)*SIN(THT))/DEN 
SUBDEN = ( .85-EB)/(l+BETA) 

EF = 1. + 2*EQ*SUBDEN 

D: = 1. +EH*SUBDEN 

EZ = .5+3.62*ET*SUBDEN 

PTOK = 24*E*Y1I/(DH**3*EL*(1-P2)) 

PV/IND = 2 .6*E’=>'(T/BH)**2.5/(EL/DM-.45*(T/DM)’V*.5) 

DENONE = (l+( .85*BETA/EB))/(1+BETA) 

DENUM = 2*T*Y/DM 

PY88 = 24*E*YBI/(DB**3*B*DEN0NE) 

• SI92 = 1/(.5+1.81*EK*SUBDEN) 

SI 92a =: 1 ./(1+EH*SUBDEN) 

SIWENK = 1./(EF**2-EF/2+.25)**.5 

PLUNCH = 2*t*Y/D/(EX**2-E>:*EZ+EZ**2)**.5 

PBRY = 2*E*T/DH/850. + 64*e*Y 1 l/(DM**3*EL) 

SIGOOR = PC •D’’'B*DEK0NE/(AB*2) 

SIGOOR = SIGOOR + E*( (D-DB)/2 . )*4*(PC/(PBRY-PC) )/R0NE**2 
SIGF = 1000*DM/(4*T) 
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SIGT = SIGF*2*(1 ,-i-2*EQ*SUBDEN) 

SIGBT = SIGF*2*ET’'SUBDEN 
SIGFOT = S IGF/S IGT 
SIGBOT = SIGBT/SIGT 
R = D/2. 

PI = 3.14159 

VT = FI*(EL*(R**2-(R-T)**2+B*((R-T)**2-(R-T-W)**2) 

1 +F* ( (R -T-W )**2 - (R-T-W-T2 )**2 ) ) 

V/T = VT-.284*12/EL 

VOL = F1*(R/12 .)**2Vfi2 

SF92A = S I 92A*DENUM*VOL/(WT* 100000) 

SFV/ENK = SIWEKK*DEKUJ1*V0L/(WT* 100000) 

Write statements for whatever solutions are required 
500 SR = SR+DSR 
700 PC = PC+100.0 

Format statements 
END 



FUNCTION SUBROUTINES 



FUNCTION YENTRD(EL,T.B,W,T2,F) 

902 TOP = 2*F*T2*W + F*T/**2 + B*W**2 - EL^'T**2 

903 BOT = 2*(F*T2 +B*W + EL*T) 

904 YENTKD = TOF/BOT 
RETURN 

END 



FUNCTION YHOM(EL,TLBLV/LT2 ,F,Y) 

922 FIRST = (EL*T**3 + B’’nv**3 + F*T2**3)/12 

923 SECOND = EL*T*(T/2+Y)**2 + B*W*(Y-W/2 )**2 

924 THIRD = F*T2*(T2/2+W-Y)**2 

925 YMOM = FIRST + SECOND + THIRD 
RETURN 

END 
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PRINTOUT OF PROGRAM 2 (Sample, continued on next page) 



T IS 3.0 INCHES AND T/DM IS 0.00529 
LAMBDA, EL, FRAME FACTOR 0.380 12.883 4 

SECOND MOMENTS :F,B,1 200.436 439.394 1481.279 

F AREA WINDENBG FORM88 FORM92 FORM92A V/ENK 

28.2750 -28193.383 590.321 1.857 1.488 1.665 

LAMBDA , EL, FRAME FACTOR 0.400 14.275 4 

SECOND MOMENTS :F,B,1 225.964 485.456 1680.225 

F AREA V/INDENBG FORM88 FORM92 FORM92A WENK 

30.022 -42473.102 609.692 1.830 1.469 1.651 

LAMBDA , EL, FRAME FACTOR 0.420 15.738 4 

SECOND MOMENTS ;F,B,1 253.112 533.708 1893.137 

F ARAE WINDENBG FORM88 FORM92 FORM92A WENK 

31.774 -90844.625 628.875 1.801 1.451 1.637 

LAMBDA , EL, FRAME FACTOR 0.440 17.273 4 

SECOND MOMENTS :F,B,1 281.881 584.134 2120.124 

F AREA WINDENBG FORM88 FORM92 FORM92A WENK 

33.531 467240.063 647.914 1.770 1.434 1.624 

U^MBDA , EL, FRAME FACTOR 0.460 18.879 4 

SECOND MOMENTS:F,B,l 312.279 636.729 2361.309 

F AREA WINDENBG FORM88 FORl'192 FORM92A WENK 

35.293 62893.883 666.852 1.736 1.420 1.612 

LAMBDA , EL, FRA MiE FACTOR 0.480 20.556 4 

SECOND MOMENTS:F,B,l 344.308 691.478 2616.743 • 

F AREA WINDENBG FOR1-I88 FORM92 FORM92A WENK 

37.058 33034.949 685.733 1.700 1.400 1.599 

LAMBDA , EL, FRAME FACTOR 0.500 22.305 4 

SECOND MOMENTS:F,B,l 377.968 748.369 2886.512 

F AREA WINDENBG FORM88 FORM92 FORM92A WENK 

38.828 22097.688 704.602 1.662 1.383 1.587 

LAMBDA , EL, FRAME FACTOR 0.520 24.125 4 

SECOND MOMENTS :F,B,1 413.264 807.395 3170.714 

F AREA WINDENBG FORM88 FORM92 FORM92A WENK 

40.600 16434.441 723.510 1.622 1.367 1.575 

LAMBDA , EL, FRAME FACTOR 0.540 26.016 4 

SECOND MOMENTS :F,B,1 450.198 868.543 3468.377 

F AREA WINDENBG FORM88 FORM92 FORM92A WENK 

42.376 12977.973 742.510 1.581 1.350 1.563 

LAl-IBtA ,EL,FRAI>1E FACTOR 0.560 27.979 4 

SECOND MOMENTS ,F,B,1 488.769 931.803 3782.55 

F AREA WINDENBG FORM88 FORM92 FORM92A WENK 

44.154 10653.945 761.666 1.538 1.333 1.551 
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FLANGE, WEB; LENGTH, THICKNESS 9.80 1.96 7.84 1.15 

LUNCHICK RATIOS :F/T,B/T .750 .0080 

LUNCHICK BRYANT OUT ROUND WT/FT EF92 EFWENK 

1.629 2478.0 54958.6 26118.07 1.246 1.327 

FLANGE, WEB: LENGTH, THICKNESS 10.10 2.02 8.08 1.19 

LUNCHICK RATIOS ,F/T,B/T .742 .0095 

LUNCHICK BRYANT OUT ROUND WT/FT EF92 EFWENK 

1.610 2526.2 55407.4 25657.38 1.253 1.339 

FLANGE, WEB: LENGTH, THICKNESS 10.39 2.08 8.31 1.22 

LUNCHICK RATIOS :F/T,B/T - .735 .0112 

LUNCHICK BRYANT OUT ROUND WT/FT EF92 EFWENK 

1.592 2571.9 55852.6 25234.33 1.257 1.348 

FLANGE, WEB: LENGTH, THICKNESS 10.68 2.13 8.54 1.26 

LUNCHICK RATIOS :F/T,B/T .727 .0131 

LUNCHICK BRYANT OUT ROUND WT/FT EF92 EFWENK 

1.574 2615.4 56292.1 24844.41 1.265 1.356 

FLANGE, WEB: LENGTH, THICKNESS 10.95 2.19 8.76 1.29 

LUNCHICK RATIOS :F/T,B/T .720 .0151 

LUNCHICK BRYANT OUT ROUND WT/FT EF92 EFWENK 

1.556 2656.7 56722.0 24483.86 1.267 1.363 

FLANGE, WEB: LENGTH, THICKNESS 11.22 2.25 8.98 1.32 

LUNCHICK RATIOS :F/T,B/T .713 .0174 

LUNCHICK BRYANT OUT ROUND WT/FT EF92 EFWENK 

1.538 2695.9 57140.1 24149.51 1.268 1.369 

FLANGE, WEB: LENGTH, THICKNESS 11.49 2.30 9.19 1.35 

LUI'ICHICK RATIOS :F/T,B/T .706 .0198 

LUNCHICK BRYANT OUT ROUND WT/FT EF92 EFWENK 

1.519 2733.4 57544.6 23838.47 1.269 1.374 

FLANGE, WEB: LENGTH, THICKNESS 11.75 2.35 9.40 1.38 

LUNCHICK RATIOS :F/T,B/T .700 .0225 

LUNCHICK BRYANT OUT ROUND V/T/FT EF92 EFWENK 

1.501 2769.0 57932.6 23548.48 1.269 1.378 

FIANGE, WEB: LENGTH, THICKNESS 12.00 2.40 9.60 1.41 

LUNCHICK RATIOS :F/T,B/T .694 .0253 

LUNCHICK BRYANT OUT ROUI'JD WT/FT EF92 EFWENK 

1.482 2803.1 58301.6 23277.38 1.267 1.379 

FLANGE ,WEB: LENGTH, THICKNESS 12.25 2.45 9.80 1.44 

LUNCHICK RATIOS :F/T,B/T .687 .0284 

LUNCHICK BRYANT OUT ROUND WT/FT EF92 EFWENK 

1.462 2835.6 58649.1 23023.39 1.265 1.381 
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START ^ 



MF, D, PC 




CALCULATE T, DM, TOD 
. l-TRITE T, TOD 



r 



CALCULATE EL, 


RATIO 




FRAMEl, DFRP, 


DFRM, NUM, 


NEND 





lA 



IS YF<^ 

NUM=NEND 



2A 



NO 



CALCULATE F, B, T2 , W, 
YB, YPI, DBT, RAT 



IS ■ 
RAT<" RATIO 



YES 



NO 

is " - 
RAT>RAT10 



NO 



FLOW CHART 
PROGRAM 2 



3 



-94- 





lA 



Loop 2 to 2A is the same as 1 to lA in its method 
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N OMENCL ATURE OF PRO GRA M 2 



B 

D 

DSR 

E 

EL 

F 

FRAMDl 

MF 

NSR 

P 

PBRY 

PC 

P LUNCH 

PTOK 

PWIND 

PY88 

R1 

R2 

R3 

SF92A 

SFV7ENK 

SIGFOT 

SIGBOT 

SI92 

SI92A 

SIV7ENK 

SR 

T 

T2 

W 

WT 

VOL 

YENTRD 

YMOM 



Width of frame in contact with shell 
Outside diameter of shell 
Incremental change in slenderness ratio 
Young's Modulus 

Length of shell between frames 
Frame overall width 
Frame overall depth 
Frame Factor 

Number of increments of slenderness ratio 
Poisson's Ratio 

Collapse pressure due to Bryant's Formula 
Required collapse pressure 
Collapse pressure due to Lunchick 
Collapse pressure due to Tokugawa , 

Collapse pressure due to Windenburg 
Collapse pressure due to' #88 
Frame v/idth/Frame depth 
Web thickness/Frame Depth 
Flange thickness/Web thickness 
Efficiency due to #92A 
Efficiency due to Wenk's Formula 
op/oT for Lunchick Formula 
obt/'^T Lunchick Formula 

^ due to #92 
tp due to #92A 
due to Wenk 
Slenderness Ratio 
Thickness of shell 
Thickness of flange 
VJidth of web 

Weight of vessel per foot length 
Volume displaced per foot length 
Subprogram to find neutral axis 
Subprogram to find I 
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APPENDIX G 
COMPUTED DATA 
and 

EXPERIMENTAL DATA FROM THE LITERATURE 
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TABLES OF MAXIMUM VALUES FROM FORMULA #92A AND MAXIMUM 



PM'iC I P LB ' S TRES s' cri terion" 



lOOh/Dj^ 



Frame Factor 4 .251 

.313 

.376 

.439 

.503 

.566 

.629 

.692 

.756 

.819 

.883 

.946 

1.010 

Frame Factor 5 .251 

.313 

.376 

.439 

.503 

.566 

.629 

.692 

.756 

Frame Factor 6 .251 

.313 

.376 

.439 

.503 

.566 

.629 

.692 

.756 

.819 

.883 

.946 

1.010 



^ max/l .2 X 



.9086 


.80 


1.785 


.9490 


.72 


1.664 


.9806 


.66 


1.574 


1.0052 


.61 


1.509 


1.0253 


.58 


1.447 


1.0426 


.54 


1.410 


1.0570 


.51 


1.375 


1.0685 


.49 


1.341 


1.0800 


.46 


1.321 


1.0901 


.45 


1.291 


1.0987 


.44 


1.266 


1.1074 


.42 


1.252 


1.1146 


.40 


1.240 


.9317 


.80 


1.716 


.9706 


.72 


1.602 


.9994 


. 6 6 


1.519 


1.0224 


.61 


1.458 


1.0411 


.57 


1.409 


1.0570 


.54 


1.367 


1.0699 


.51 


1.334 


1.0814 


.49 


1.304 


1.0916 


.47 


1.278 


.9504 


.82 


1.642 


.9864 


.72 


1.554 


1.0138 


.65 


1.485 


1.0354 


.61 


1.419 


1.0541 


.57 


1.373 


1.0685 


.54 


1.334 


1.0814 


.51 


1.305 


1.0916 


.49 


1.276 


1.1016 


.47 


1.253 


1.1102 


.45 


1.234 


1.1175 


.43 


1.219 


1.1246 


.42 


1.201 


1.1304 


.40 


1.191 
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TABLES O F M AXIflUM VALUES FROM WENKS FOR.MULA 
HENCKY-yp_N_~MISES YIELD CRITERION-r.IDPLANE , MIDBAY 



100h/Dj„ 



Frame Factor 4 .376 

.439 
. 503 
.566 
.629 
.692 
.756 
.819 
.883 
.946 
1.010 

Frame Factor 5 .376 

.439 

.503 

.566 

.629 

.692 

.756 

Frame Factor 6 .376 

.439 

.503 

.566 

.629 

.692 

.756 

.819 

.883 

.946 

1.010 



Tj max/1 ,2 


A 




1.1333 


.82 


1.646 


1.164 


.76 


1.611 


1.198 


.72 


1.562 


1.214 


.67 


1.533 


1.237 


.63 


1.506 


1.248 


.60 


1.481 


1.265 


.58 


1.454 


1.270 


.55 


1.438 


1.280 


.52 


1.425 


1.292 


. 50 


1.409 


1.303 


.49 


1.391 


1.1592 


.84 


1.620 


1.1880 


.76 


1.571 


1.2139 


.70 


1.538 


1.2341 


.67 


1.497 


1.2500 


.63 


1.472 


1.2644 


.60 


1.448 


1.2759 


. 57 


1.428 


1.1794 


.84 


1.573 


1.2068 


.76 


1.540 


1.2298 


.71 


1.502 


1.2471 


.66 


1.476 


1.2609 


.62 


1.451 


1.2744 


.60 


1.423 


1.2859 


.57 


1.405 


1.2960 


.54 


1.390 


1.3032 


.52 


1.374 


1.3118 


.51 


1.3 56 


1.3176 


.49 


1.345 



FPIAME PARA-METERS AT MAXIMUM EFFICIENCY #92A 



lOOh/D^ 

Frame Factor 4 .251 

.313 

.376 

.439 

.503 

.566 

.629 

.692 

.756 

.819 

.883 

.946 

1.010 

Frame Factor 5 .251 

.313 

.376 

.439 

.503 

.566 

.629 

.692 

.756 

Frame Factor 6 .251 

.313 

.376 

.439 

.503 

.566 

.629 

.692 

.756 

.819 

.883 

.946 

1.010 



Il/Ip 


lOOb/D 


4.758 


.2195 


5.335 


.2350 


5.860 


.2485 


6.340 


.2640 


6.817 


.2720 


7.238 


.2780 


7.655 


.2850 


8.072 


.2930 


8.464 


.2960 


8.871 


.3045 


9.267 


.3120 


9.662 


.3145 


10.065 


.3150 


5.064 


.2065 


5.672 


.2215 


6.224 


.2340 


6.735 


.2420 


7.212 


.2525 


7.679 


.2610 


8.124 


.2675 


8.567 


.2745 


9.001 


.2810 


5.364 


.1990 


5.958 


.2110 


6.533 


.2220 


7.062 


.2315 


7.566 


.2400 


8.056 


.2475 


8.527 


.2535 


8.995 


.2605 


9.451 


.2660 


9.919 


.2700 


10.391 


.2725 


10.853 


.2785 


11.354 


.2790 



80 

72 

66 

61 

58 

54 

51 

49 

46 

45 

44 

42 

40 

80 

72 

66 

61 

57 

54 

51 

49 

47 

82 

72 

66 

61 

57 

54 

51 

49 

47 

45 

43 

42 

40 
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FFiAME PARAMETERS AT MAXIMUM EFFICIENCY #WENK 



lOOh/D^ 



Frame Factor 4 .376 

.439 

.503 

.566 

.629 

.692 

.756 

.819 

.883 

.946 

1.010 

Frame Factor 5 .376 

.439 

.503 

.566 

.629 

.692 

.756 

Frame Factor 6 .376 

.439 

.503 

.566 

.629 

.692 

.756 

.819 

.883 

.946 

1.010 





Ii/If 


lOOb/D 


82 . 


6.192 


.2800 


76 


6.647 


.2925 • 


72 


7.085 


.3055 


67 


7.468 


.3140 


63 


7.843 


.3210 


60 


8.213 


.3290 


58 


8.580 


.3375 


55 


8.930 


.3420 


52 


9.280 


.3440 


50- 


9.634 


.3495 


49 


9.986 


.3555 


84 


6.568 


.2670 


76 


7.008 


.2755 


70 


7.433 


.2830 


67 


7.878 


.2950 


63 


8.264 


.3010 


60 


8.657 


.3090 


57 


9.041 


.3140 


84 


6.853 


.2540 


76 


7.311 


.2620 


71 


7.768 


.2715 


66 


8.195 


.2780 


62 


8.616 


.2835 


60 


9.037 


.2930 


57 


9.446 


.2980 


54 


9.859 


.3010 


52 


10.273 


.3060 


51 


10.678 


.3130 


49 


11.106 


.3160 
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EXPERIMENTAL 


DATA FROM 


REFERENCE 


9 






Model Number 


lOOh/D 


^Jy/1000 


X 


V^92A 


V^exp 


1 


.88 


46.8 


.41 


1.43 


1.68 


2 


.784 


47.5 


.42 


1.42 


1.74 


BR7m 


.784 


59.2 


.46 


1.44 


1.62 


BR7 


.599 


44.8 


.59 


1.11 


1.3 


5 


.814 


47.1 


.61 


1.01 


1.24 


6 


.474 


49.6 


.67 


1.09 


1.24 


7 


. 581 


47.0 


.70 


.997 


1.11 



Y'' exp ^represents y/ determined using actual collapse 
pressure and determined using collapse pressure 

calculated by formula #92A. 
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EXPERIMENTAL DATA FROM REFERENCE 19 



General Model Parameters 
Left/D^j = .1405 
E = 30,000,000 Ib/in^ 



Light ] 
MODEL 1 


Frames 

NUMBER 


10 Oh/D 


CTy/1000 


X 


V'eKp 


79 




.301 


27 


.876 


.92 


82 




.303 


28 


.855 


.90 


76 




.307 


42 


1.06 


.66 


88 




.293 


27 


.890 


.85 


75 




.294 


39 


1.06 


.91 


81 




.295 


28 


.900 


.93 


84 




.300 


29 


.909 


.81 


87 




.300 


29 


.905 


.81 


80 




.310 


28 


.840 


.87 


85 




.311 


29 


.884 


.83 


92 




.288 


32 


1.00 


.85 


91 




.290 


32 


.975 


.84 


93 




.308 


29 


.886 


.88 


96 




.298 


28 


.896 


.85 


94 




.298 


28 


.896 


.84 


90 




.300 


27 


.88 


.82 


In the 


computation of ' 


\ values , 


h/D was used 


instead 


of h/D^ resulting in a 


small error which was 


not criti( 



in the graphical methods used. 



represents the ratio of actual collapse pressure to 



the hoop stress collapse pressure 
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ALUM I NAUT 



In the following section, A, ip, and n have 
been found for the pressure hull in Aluminaut. The 
dimensions of that craft were taken from References 
(5) and (15) and are for the cylindrical hull alone. 



Operating Depth 
Safety Factor 



E 

P 



15,000 ft. 


LOA = 


33.3 


ft 


1.4 


tl 

Q) 

1^ 


6.6 


ft 


60,000 psi 


D 


00 

o 


ft 


29,000,000 psi 


Dm = 


7.5 


ft 


.1 Ib/in^ 


h 


.5 


ft 




2 . 






Pc 

oh 

1.25 



p^ = 1.4 X 15,000 X .446 
P(^ = 9400 psi 



3 . n 



Pc^ 

W 10“ 



n = 2.65 in. 



From the above computations, it is apparent that the 
designers of Aluminaut produced a vessel whose parameters 
put it in the low lambda design region, Hov/ever, it has 
been noted (15) that a more efficient structure is possible. 



I 
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